Maximally discordant mixed states of two qubits 
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We study the relative strength of classical and quantum correlations, as measured by discord, for 
two-qubit states. Quantum correlations appear only in the presence of classical correlations, while 
the reverse is not always true. We identify the family of states that maximize the discord for a given 
value of the classical correlations and show that the largest attainable discord for mixed states is 
greater than for pure states. The difference between discord and entanglement is emphasized by 
the remarkable fact that these states do not maximize entanglement and are, in some cases, even 
separable. Finally, by random generation of density matrices uniformly distributed over the whole 
Hilbert space, we quantify the frequency of the appearance of quantum and classical correlations 
for different ranks. 

PACS numbers: 03.67.Mn, 03.65.Ta 



I. INTRODUCTION 

One of the most striking features of quantum mechan- 
ics is entanglement, first considered (although not by 
that name) by Einstein, Podolsky, and Rosen in their 
seminal paper in 1935 [1 . This is an exclusively quan- 
tum feature of composite states that can not be written 
as mixtures of product states. Theoretical and experi- 
mental research activity to characterize entanglement has 
been particularly intense in the last decade (see review [2] 
and references therein) , being part of a broader endeavor 
to explore distinctive aspects of quantum versus classi- 
cal physics and novel resources for quantum information 
purposes [3 . An important issue considered by several 
authors [4-6 is the existence of quantum correlations be- 
yond entanglement, in separable states. As a matter of 
fact, examples of improved quantum computing tasks not 
relying on entanglement have been reported [7^. 



II. QUANTUM CORRELATIONS: THE 
DISCORD 



Two complementary approaches on quantum correla- 
tions are receiving great attention [4, 5 . In Ref. i4|, 
quantum correlations (quantum discord) have been as- 
sociated to the difference of two classically equivalent 
expressions for the mutual information, X and J . In 
particular, the quantum mutual information is defined 
as X{q) = S{qa) + S{qb) — S{g)^ where S stands for 
the von Neumann entropy and Qa^b) is the reduced den- 
sity matrix of each subsystem. The classically equiva- 
lent expression stemming from Bayes rule is J{g){u^y = 

S{qa) — S{A\{II^}), with the conditional entropy defined 
as S{A\{Uf}) = E^P^S{QA\uf). Pi = TiABiUfg), and 
where ^^in^ = ^f Q^f /Pi is the density matrix after a 
complete projective measurement ({11^}) has been per- 
formed on B. Quantum discord is obtained minimizing 



the difference X{g) — J{q)'. 

Sa:b{q) = min [SiQB) - S{g) + S{A\{Uf})] , (1) 
{nf } 

that is, when measurement is performed in the basis 
which disturbs the state the least. A complementary ap- 
proach was described in Ref. [5j , defining classical corre- 
lations and showing that total correlations given by the 
mutual information are actually larger thAn the sum of 
the classical correlations and entanglement E [8 . As a 
matter of fact, the quantum mutual information can be 
seen as the sum of quantum correlations [4 Sa:b{q) and 
the classical correlations [5 maxr^Bi J'(^)rnB|. We re- 
mark that, in [4 , the discord is defined in terms of orthog- 
onal (perfect) measurements. Even if possible generaliza- 
tions to positive-operator- valued measurements (POVM) 
were considered at the end of that paper, as well as in 
^ , calculations of discord in the literature generally con- 
sider only orthogonal measurements (see, e.g., [4[ 13 |9]- 

In contrast with state separability, this new paradigm 
of quantumness of correlations is measurement oriented, 
considering an experiment where all information of a 
system A is extracted by measuring another system B. 
According to this measure, a state is classically corre- 
lated only when consecutive measurements of system B 
yield the same picture of the state of system A, which 
is achieved after decoherence into the pointer basis of 
B [H [9]. For pure states, quantum discord is equiva- 
lent to entanglement and actually has the same value of 
classical correlations [5 . On the other hand, when mixed 
states are considered, entanglement does significantly de- 
part from the quantum discord, the difference being pos- 
itive for some states and negative for others [TO] . As the 
definition of discord comes from a minimization over all 
possible measurement basis, only a few general results 
have been reported. Analytic expressions are known for 
states of two qubits with maximally mixed marginals [10], 
for X-shaped states [H , and also for Gaussian states of 
continuous variable systems [12 . 



In this paper, we explore the whole Hilbert space of 
two qubits to gain insight on their correlations for mixed 
states of different ranks. Our main goal is to discern 
the proportion of quantum to classical correlations be- 
tween the two qubits. We find the most nonclassical two- 
qubit states, i.e., the family with maximal quantum dis- 
cord versus classical correlations, were formed by mixed 
states of rank 2 and 3, which we name maximally dis- 
cordant mixed states (MDMS). The analogous effort to 
identify largest deviations from classical states has led 
to mixed states maximizing entanglement versus purity 
[13 . In contrast with maximally entangled mixed states 
(MEMS), where a geometrical property such as separa- 
bility could be considered with several constraints, the 
MDMS are naturally defined and allow us to quantify 
the relative strength of quantum and classical correla- 
tions, which are related in a closed from. The MDMS 
proposed here do not maximize entanglement for a given 
amount of classical correlations; part of them are, in fact, 
separable. This pinpoints the fundamental difference be- 
tween entanglement and discord for mixed states, in op- 
position to their exact equivalence for pure states. Their 
discord also implies that quantum correlations are always 
accompanied by classical correlations, while the reverse 
is not always true. Furthermore, we study the proba- 
bility of states with a given amount of discord in the 
whole two-qubit Hilbert space, supporting the recent re- 
sult that the closed set of purely classically correlated 
states {6a:B = 0) has measure zero [14 , and we com- 
pare it with the probability for classical correlations and 
entanglement. Only such states with no discord have 
been shown to ensure a future non-negative evolution in 
the presence of dissipation [15] , while discord can not be 
made zero in a finite time by any Markovian map [M]. 
Furthermore, even in the presence of a noisy environ- 
ment, for some family of initial states, discord can be 
robust under decoherence for a finite time [16 . Experi- 
mental results with polarization entangled photons also 
have been reported recently [ 17] . 



III. MIXED STATES WITH LARGEST 
DISCORD 

As pure states with maximum entanglement that are 
also maximally discordant. Bell states are a natural start- 
ing point to identify states with large discord; thus, by 
mixing these states with other components, it is quite 
plausible that we find states with a large proportion of 
quantum versus classical correlations. We then consider 
an example of a mixture of any Bell state lip) with an- 
other orthogonal pure state, i.e., g = e|'0)('0| + (1 — 
e)|0)(0|. If \(j)) is any other Bell state, then Sa-.b = E 
(with E being entanglement) and J7 = 1, and a worse 
discord is found than in that of pure states with the same 
classical correlation J. In contrast, a mixture of a Bell 
state and a state of the computational basis of the op- 
posite parity sector, gives a huge amount of discord. We 



thus consider states 



^ = e|$+)($+| + (l-e)|01)(01| 



(2) 



with !<!>+) = (|00) + |11))/V2 as usual. With local uni- 
tary operations, which leave discord invariant, we can 
obtain from this ansatz any combination of a Bell state 
mixed with a computational basis state of opposite parity 
(number of I's in the state). The expression of discord 
for states ^ is invariant under permutation of the indi- 
vidual labels A ^ B. As a matter of fact, we find that 
states Q maximize the symmetrized version of discord 
{^A:B + ^B:a)/2 for all rank-2 matrices. 
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FIG. 1: (Color online) Quantum discord (Sa-.b) versus clas- 
sical correlations (J') for two-qubit states. The MDMS fam- 
ily (continuous line) gives two segments for rank 3 (g^^^-') 
and one for rank 2 (g^^'^-'). Layers of 10^ random matrices 
of rank 2 (dark points), 3 (intermediate color), and 4 (lighter 
color) are superimposed. For pure states (dotted-dashed line), 

Sa:B =J = E. 



It can be shown that, when we compare with a numer- 
ical scan of Hubert's space, this state is too symmetric. 
In fact, a better option in terms of discord is obtained if 
some amount of entanglement is sacrificed for the good 
of quantum correlations. This results from asymmetriz- 
ing the maximally entangled (Bell) state leading to the 
ansatz 



,(R2) 



e|<l+)(l>+| + (l-e)|01)(01|. 



(3) 



with |<l+) = v^|00) 
a Bell state for p 



vr 



p|ll)), which coincides with 
1/2. The increase of discord for 
the states ([s]) with respect to Q highlights the im- 
portance of the asymmetric definition of quantum dis- 
cord, based on the asymmetric operation of measuring 
B in order to know about A. The discord for this fam- 
ily can be written once we know the conditional en- 
tropies min|nf}5'(A|{nf}) = (xlog2 j^ -log2^)/2, 
with X = y^l — Ay and y = e(l — p){l — e), while for 
5b -.A we need to use y = ep(l — e). The total and reduced 



entropies are easy to calculate once we notice that the 
ansatz is given in spectral decomposition, though we do 
not give the whole expression for reasons of space. 

The family of MDMS is obtained for an optimal func- 
tion eopt.(p) through the use of Lagrange multipliers, as 
detailed later in this paper. Once this optimal curve 
eopt.(p) is used, Eq. ([3| gives the family of states that 
maximize the quantum part of correlations for a given 
classical part , when all rank 2 (R2) states are consid- 
ered. As shown in Fig. [ij these ^^^^^ states are the 
MDMS for a large range of classical correlations. 

In order to find the states that maximize Sa-.b^ we 
also need to use rank-3 states. In this case, it can be 
checked that asymmetrization of the Bell-state compo- 
nent |(/)) does not help, and that the best choice for an 
ansatz is 
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^(R3) ^ e|$+){$+| + (l-e)(m|01)(01| + (l-m)|10)(10|). 

(4) 
As before, any combination of Bell state plus two compo- 
nents of opposite parity belonging to the computational 
basis, will do. The optimal eopt.(^Ti) is discussed later 
and leads to the family of states ^^^^^ maximizing dis- 
cord for small classical correlations. This optimal fam- 
ily has the property of being separable (not entangled), 
as shown in Fig. [2J while it maximizes quantum dis- 
cord, highlighting the inequivalence of these measures of 
quantumness. It is actually found that, for these states, 
the discord amounts to the weight of the Bell compo- 
nent, and the simple relation Sa-.b = ^b-.a = ^ holds. 
For completeness, the entanglement (as quantified by 
the concurrence [8 ) yields E{g^^'^'^) = 2ey/p{l — p) and 
E{gi^^)) = max[0,e - 2(1 - e)v'm(l-m)]. Although 
MDMS of rank 3 are separable, the MDMS of rank 2 have 
a high amount of entanglement, even if they do not max- 
imize it (Fig. [2]). As mentioned before, asymmetrization 
of the Bell component increases quantum correlations at 
the expense of entanglement. 

An intriguing state is that of the singular point for 
^^^3) shown in Fig. [ij ^^^^p = (|<|>+)(<|>+| + |01)(01| + 
|10)(10|)/3 reaches the lowest possible purity for a rank- 
3 state and is separable, yet with a high level of discord. 
Another important feature emerging from Fig. [l] is that 
MDMS have a discord larger than pure states {Sa-.b > sJ) 
and satisfy J = only when Sa-.b = 0, thus showing the 
lack of states with finite quantum without classical corre- 
lations. In other words, no state of two qubits is purely 
quantum. Maybe less surprisingly, there are no states 
with finite entanglement and zero classical correlations 
(Fig.i. 



IV. DETAILED ANALYSIS 

We first consider the commonly accepted definition 
of discord, obtained expressing the measurement pro- 
jectors as Uf = \^f){^fl j = 1,2, with |V^f) = 

-e-^^sin(9|0)+cosl9|l) 



FIG. 2: (Color online) Entanglement (E) vs classical correla- 
tions (sJ). Random matrices and lines for the MDMS family 
as in Fig. [ij MDMS of rank 3 are separable, while g^^'^-' states 
show large entanglement, even if not maximum, as seen in the 
inset. 



[4 . The angles minimizing Eq. (11) for the states 
and thus giving the correct discord, are = 7r/4 + utt/ 
and any value of (j). We then can choose (j) = and 
hence the optimal projectors are Ilf = |+^)(+^| and 

nf = |-^>(-^|,with|±) = (|o)±|i))/v^. 

The method of Lagrange multipliers allows us to maxi- 
mize first the discord for the rank-2 family Sa-.b{q''^'^^) = 
S{e^p) while keeping its classical correlations J'{g^^'^'^) = 
j7o(e,p) constant (and the same for the rank-3 family 
^(R3)j^ Xhis is achieved through definition of the func- 
tion A(e,p, A) = S{e,p) + A(j7(e,p) — jTo) where A is the 
Lagrange multiplier, and J'o is an arbitrary but fixed 
amount of classical correlation. The extremization pro- 
cedure is then simply the simultaneous solution of the 
three equations 9^ A = 0, with /i = A, e,p. From the first 
equation, J{e^p) = Jq is obtained, as expected. From 
the last two equations, we can isolate A yielding the ex- 
tremality condition 



d,S/d,J = dpS/dpJ. 



(5) 



cosi9|0) 



I^ - 
' sin^ 



|1> and \^i) 



We stress that this condition is equivalent to maximiza- 
tion of S versus X, or minimization of J with respect to 
X, due to the closed relation X = S -\- J'. These quanti- 
ties present nontrivial trigonometric relations, leading to 
a transcendental equation the solution of which can only 
be given numerically. 

The same procedure is followed for the rank-3 family 
^(R3)^ with m playing the role of p. In this case, obtain- 
ing 6opt.(^) is a bit trickier, due to the fact that there 
are two optimal angles (each of them good for different 
ranges of e and m)^ = 0,7r/4; the angle (j) is again 
non important. We can consider for the moment that 
projector maximization of discord has been simplified to 
^A:b(^^^^^) = ^^{So^S^/4^)^ the latter being functions of 
e and m. The goal is to find the zero(es) of the function 



d^5/deJ — dm^/dmJ^ which of course needs the knowl- 
edge of when to use one angle or the other. However, the 
problem is greatly reduced by noticing that the latter 
function is positive when using 5q and negative when us- 
ing ^7^/4- This means that the zero of such function occurs 
exactly (and conveniently) when (5o(e,m) = (57^/4(e,m). 
Again, the solution to this transcendental equation can 
only be given numerically. 

Finally, the MDMS are a family of states 

^MDMS ^ e|l>+)(l>+| + (l-e)(m|01)(01| + (l-m)|10)(10|), 

(6) 
where the optimum choice of parameters gives the three 
curves in Fig. [l] two of them rank 3 and the other rank 
2. The first curve, going from zero discord up to the 
cusp, is the rank 3 family g^^^^ with eopt.(^) given by 
the solution of 5q = S^/^. It is restricted to the domain 
m G [0, 1], e G [0, 1/3]. The second branch of MDMS is 
given by ^^^^^ with m = 1/2 with domain e G [1/3, 0.385], 
approximately. These two curves correspond to separa- 
ble states, as noted above (Fig. [2]). The remaining curve 
of MDMS is the rank 2 family g'^^'^^ when the optimal 
function eopt.(p) given by Eq. (I5| is used, and for e ap- 
proximately in the interval [0.408, 1]. One might wonder 
how the picture changes if more general (nonorthogonal) 
measurements are considered. It has been shown that, 
for two qubits, discord is extremized exclusively by rank- 
1 POVMs with a maximum of four elements [211 l22] . 
Perfect orthogonal measurements correspond to the case 
with the two elements considered above. Considering 
POVMs with measurement operators Ei^ the measured 
density matrix takes the form gA\Ei = ^iQ/Pi with prob- 
ability Pi = Ti AB{Eig). Even using the general measure- 
ment given by POVMs of four elements, we find the same 
discord for the MDMS, meaning that they represent the 
absolute border of maximally nonclassically correlated 
states of two qubits. A detailed analysis about the full 
Hilbert space will be presented elsewhere [23^ . 



to lower ranks, we observe that rank 2 matrices yield 
~ 10.76%, and rank 3 yield ~ 16.3%. Finally, as shown 
in the insets in Fig. [3J the border of MDMS seems to be 
rather improbable to find in the space of two qubit states, 
except for the middle branch in the cusp, meaning that 
such extremely nonclassical states are quite rare. 



pisy 




FIG. 3: (Color online) Probability (density) to find a two- 
qubit state with a given amount of quantum discord 6, classi- 
cal correlations JT, and entanglement E respectively, for ranks 
2 (dashed), 3 (dot-dashed), and 4 (dot-dot-dashed). The in- 
sets show these probabilities (larger for light color) for the 
quantity under study {y axis) against classical correlations J' 
{x axis), for different ranks. 



VI. CONCLUSIONS 



V. STATISTICS 

Since our random generation of density matrices is uni- 
form in the Hilbert space, preserving the Haar measure 
[IHl [19] , we can measure the frequency of the appearance 
of states with different properties, as shown in Fig. [3] for 
different ranks. Some main features arise for all quan- 
tities investigated: (i) Zero correlations (be they quan- 
tum, classical, or mutual information) have zero probabil- 
ity. Notably, entanglement is the exception (consistently 
with Ref. [E]), where only rank-2 states have a probabil- 
ity zero of being separable, (ii) The lower the rank, the 
higher the typical amount of correlations. This is quite 
understandable, since higher ranks describe more mixed 
states, (iii) It is more probable to find states with more 
abundance of classical rather than quantum correlations. 
We note that only ~ 7.45% of the two-qubit states has 
greater discord than classical correlations. If restricted 



The unique family of two-qubit mixed states with max- 
imal proportion of quantum discord versus classical cor- 
relations, the MDMS in Eq. (IgI), has been identified. Part 
of them have rank 2 and are highly, although not max- 
imally, entangled, while the other part has rank 3 and 
is separable, thus providing another evidence of the in- 
equivalence of these two measures of quantumness. We 
have shown that the presence of discord is a sufficient 
but not necessary condition to have nonvanishing classi- 
cal correlations. The uniform generation of states (ran- 
dom states preserving Haar measure) allowed us to find 
the probabilities and typical values of classical and quan- 
tum correlations, as well as entanglement. We verified 
that completely (either quantum or classical) uncorre- 
lated states are very rare, as well as extreme nonclassical 
states. The identification of MDMS, together with the 
ability to experimentally generate [20] and characterize 
[17 these states, is a key tool to establish the fundamen- 
tal difference in performance of quantum versus classical 



information |7]. 
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